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In this short review on stellar convection dynamics I address the following, currently very topical, issues: (1) the surface 
effects of the Reynolds stresses and nonadiabaticity on solar-like pulsation frequencies, and (2) oscillation mode lifetimes 
of stochastically excited oscillations in red giants computed with different time-dependent convection formulations. 
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1 Introduction 
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Seismic data of solar-type oscillations in distant stars are 
being provided with high precision by the space missions 
CoRoT (Baglin et al. 2009) and Kepler (Borucki et al. 2009; 
Christensen-Dalsgaard 2007), and hopefully also soon from 
the ground-based Danish SONG project (Grundahl et al. 
2010). In contrast to solar observations the observed oscilla- 
tion modes in distant stars are of only low spherical degree 
I. However, these are the modes which provide information 
from the deepest layers of the star. Some theoretically im- 
portant properties, such as the gross structure of the energy- 
generating core and the extent to which it is convective, and 
the large-scale variation of the angular velocity, will become 
available. Such information will be of crucial importance 
for checking, and then calibrating, the theory of the struc- 
ture and evolution of stars. 

Solar-type oscillations are characterized by being in- 
trinsically damped and stochastically driven by the near- 
surface turbulent velocity field (e.g., Goldreich & Keeley 
1977, Balmforth 1992b, Houdek et al. 1999, Samadi et al. 
2008). The turbulent convection produces acoustic noise 
in a broad frequency range, thereby exciting many of the 
global solar-type oscillations to observable amplitudes. The 
resulting rich acoustic power spectrum makes these modes 
particularly interesting for studying fluid-dynamical aspects 
of the stellar interior, and it is our hope that their mode 
properties can be used to probe the still ill-understood near- 
surface regions. The physics of these near-surface regions 
is, besides many other effects, further complicated by the 
dynamics of the transport of the turbulent fluxes of heat 
and momentum and their coupling to the pulsations. Several 
attempts have been made in the past to model the effects 
of convection dynamics on the oscillation frequencies, for 
which a time-dependent convection formalism is required. 
The influence of the momentum flux and its pulsationally 
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induced perturbation upon adiabatic and nonadiabatic pul- 
sation frequencies were investigated first by Gough (1984), 
using a local time-dependent mixing-length formulation for 
convection (Gough 1977a), but a simplified analytical ap- 
proximation to the eigenfunctions in the atmosphere. Gough 
concluded that the inclusion of the momentum flux in the 
mean model reduces the frequency residuals between ob- 
servations and model computations, whereas the additional 
dynamical effects through their pulsationally induced per- 
turbations actually aggravates this frequency discrepancies. 
Balmforth (1992a) studied these effects in a more consistent 
way and by means of the more sophisticated nonlocal time- 
dependent mixing-length formulation by Gough (1977b). 
Balmforth also concluded that the modification of the mean 
stratification by the momentum flux substantially decreases 
the adiabatic frequency residuals, whereas the effects of the 
momentum flux perturbation and nonadiabaticity lead to an 
increase in the modelled oscillation frequencies. 

In Section 2 we discuss the results of more recent at- 
tempts to model the effects of the near-surface layers on 
the oscillations frequencies in the Sun and other stars. The 
other topical issue on the comparison between observed p- 
mode lifetimes in red giants and theoretical predictions, ob- 
tained with different time-dependent convection models, is 
addressed in Section 3. 

2 The effect of convection dynamics on the 
oscillation frequencies 

The near-surface layers in the Sun and other stars are still 
poorly understood, mainly because of the rather difficult 
to model complex physics of the interaction of convection, 
radiation, magnetism and rotation. Typical stellar structure 
calculations treat these superficial layers with simplified at- 
mosphere models and the stratification of the superadia- 
batic region by means of a mixing-length approach (e.g. 
Bohm-Vitense 1958). Furthermore, the effects of the mean 
Reynolds stresses on the hydrostatic equilibrium are essen- 
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Fig. 1 Scaled differences between observed GONG fre- 
quencies and adiabatically computed frequencies of the 
'standard' solar Model S. The degree I of the oscillation 
modes is indicated by the colour bar (from Christensen- 
Dalsgaard et al. 1996). 



tially always ignored. Linear pulsation calculations typi- 
cally adopt the adiabatic approximation and ignore the mo- 
mentum flux perturbations (turbulent pressure fluctuations). 
However, nonadiabatic effects and the fluctuations of the 
turbulent fluxes (heat and momentum) do modify the mod- 
elled pulsation eigenfunctions and consequently also the os- 
cillation frequencies (e.g., Gough 1980, Balmforth 1992b, 
Houdek 1996). 

These effects are predominantly confined to the upper-most 
stellar layers, where the vertical scale of low- or intermedi- 
ate degree modes is much less than the horizontal scale and 
when I is low the mode inertia is quite insensitive to degree 
I. Consequently the observed oscillation properties of low- 
degree modes depend little on the value of I (e.g., Libbrecht 
1986). Moreover, the influence of the near-surface effects 
on the oscillation frequencies of low-degree modes is pre- 
dominantly a function of frequency alone (e.g., Christensen- 
Dalsgaard & Thompson 1997), whereas for modes with in- 
termediate degree a weak degree dependence is observed. 
This degree dependence, however, can be described by a 
simple scaling Q„i = I n i/Io with mode inertia I n i (e.g., 
Aerts, Christensen-Dalsgaard & Kurtz 2010, §7.1.4), where 
n is the radial mode order, and Iq{v) is obtained from inter- 
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Fig. 2 Scaled adiabatic frequency differences between a 
model for which the near-surface layers were represented 
by a hydrodynamical simulation, and the 'standard' solar 
Model S (from Rosenthal et al. 1995). 



polating the radial mode inertia I„q to the nonradial oscil- 
lation frequency v n i. For a 'standard' solar model, such as 
Model S (Christensen-Dalsgaard et al. 1996), the (scaled) 
differences between the solar and model frequencies are in- 
deed dominated by the near-surface effects and are predom- 
inantly a function of frequency alone. Differences between 
frequencies observed with the GONG instruments and adi- 
abatically computed Model S frequencies are illustrated in 
Fig-D The increase of the frequency residuals with oscilla- 
tion frequency depends on the modelling details of the func- 
tional form of the acoustic cutoff frequency v &c (as it affects 
the acoustic potential) with radius in the near-surface lay- 
ers. For an isothermal atmosphere v &c = c/AttH, where c is 
the adiabatic sound speed and H the pressure scale height. 
In the Sun v ac ~ 5.5 mHz. Thus, we see from Fig. Q] that 
the influence of the near-surface layers is most important, 
when the oscillation frequency v is comparable to the acous- 
tic cutoff frequency v ac . The cutoff frequency determines 
the location at which an incident acoustic wave is reflected 
back into the stellar interior, and the lower the frequency v, 
the deeper the location at which this reflection takes place. 
For modes with frequencies v much less than ^ ac reflection 
takes place so deep in the star that the modes are essentially 
unaffected by the near-surface structure. When v is compa- 
rable with v &c , however, the inertia of the near-surface layers 
is a considerable fraction of the total mass above the reflect- 
ing layer, leading to a greater modification to the phase shift 
in the spatial oscillation eigenfunctions, and, through the 
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dispersion relation, also to a change in frequency. The in- 
ertia of the essentially hydrostatically moving near-surface 
layers depends on mass and consequently on the equilib- 
rium pressure near the photosphere. 

2.1 The effect of the Reynolds stresses in the mean 
structure 

From the discussion before we conclude that the details of 
modelling the hydrostatic equilibrium structure in the near- 
surface layers play an important role in describing the resid- 
uals between observed and modelled oscillation frequen- 
cies, particularly for modes with v close to v ac . Almost all 
stellar model calculations consider only the gradient of the 
gas pressure p g in the equation of hydrostatic support. In the 
convectively unstable surface layers, however, the turbulent 
velocity field u contributes to the hydrostatic support via 
the Reynolds stresses (puu) (angular brackets denote an 
ensemble average), the (r,r) component p t := (pu r u r ) of 
which, acts as a pressure term additionally to the gas pres- 
sure p g (e.g., Gough 1977ab). Solar models indicate that the 
turbulent pressure p t can be as large as 15% of the total pres- 
sure p = p g + pt, as illustrated in Fig.[5]for three different 
simulations and models of the Sun. 

Hydrodynamical simulations of stellar convection have 
enabled us to estimate the turbulent pressure (see Fig. [3]). 
Rosenthal et al. (1995) investigated the effect on adiabatic 
eigenfrequencies of the contribution that the turbulent pres- 
sure makes to the mean hydrostatic stratification: he ex- 
amined a hydrodynamical simulation by Stein & Nordlund 
(1991) of the outer 2% by radius of the Sun, matched con- 
tinuously in sound speed to a model envelope calculated, as 
in a 'standard' solar model, with a local mixing-length for- 
mulation. The resulting frequency shifts of adiabatic oscilla- 
tions between the simulations and the 'standard' solar refer- 
ence model, Model S, are illustrated in Fig. [2] The frequency 
residuals behave similarly to the data in Fig.Q] though with 
larger shifts at higher oscillation frequencies. 

If turbulent pressure is considered in the mean structure, 
additional assumption have to be made about the turbulent 
pressure perturbation in the adiabatic pulsation equations. 
The first adiabatic exponent 71 = (d\np g /dlnp) s (s being 
the specific entropy), is a purely thermodynamic quantity 
and is expressed by means of the gas pressure p g . Conse- 
quently, in the presence of turbulent pressure p t , such that 
the total pressure p = p g + p t satisfies the equation of hy- 
drostatic equilibrium, 71 experiences a modification of the 
form 
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and the Lagrangian perturbation in the total pressure is 

Sp = Sp E + S_p 1= ^S_p (2) 

p p p p 

Nonadiabatic pulsation calculations with the inclusion of 
the turbulent pressure fluctuation 5p t (Houdek 1996), and 
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Fig. 3 Turbulent pressure p t = (p"u%) over total pressure 
p = p g + pt as a function of the depth variable R Q — r 
for various solar simulations and models (R Q is the pho- 
tospheric solar radius). Results are shown for the nonlocal 
mixing-length model by Gough (1977ab, solid curve) and 
from hydrodynamical simulations by Trampedach (1999, 
personal communication, dashed curve) and Ludwig (2005, 
personal communication, dot-dashed curve). 



hydrodynamical simulation results (Rosenthal et al. 1995) 
indicate that Spt varies approximately in quadrature with 
the other terms in the linearized momentum equation, and 
hence contributes predominantly to the imaginary part of 
the frequency shift, i.e. to the linear damping rate. The La- 
grangian perturbation 5p g , however, responds adiabatically. 
Therefore 6pt can be neglected in equation (0, i.e in the 
calculation of the real adiabatic eigenfrequencies it is as- 
sumed that Sp/p ~ Spg/p ~ jiSp/p. With this assumption 
71 ~ (p g /p)ji, and the only modification to the adiabatic 
oscillation equations is the replacement of 71 by 71 (Rosen- 
thal et al. 1995). 

2.2 The effects of nonadiabaticity and momentum flux 
perturbation 

The effects of nonadiabaticity and convection dynamics on 
the pulsation frequencies were, for example, studied by 
Balmforth (1992b), Rosenthal et al. (1995) and Houdek 
(1996). In these studies the nonlocal, time-dependent gen- 
eralization of the mixing-length formulation by Gough 
( 1977ab) was adopted to model the heat and momentum flux 
consistently in both the equilibrium envelope model and in 
the nonadiabatic stability analysis. Houdek (1996) consid- 
ered the following models: 

L.a A local mixing-length formulation without tur- 
bulent pressure pt was used to construct the 
mean envelope model. Frequencies were com- 
puted in the adiabatic approximation assuming 
Spt = (see Section l2~TI ). 
NL.a Gough's (1977ab) nonlocal, mixing-length 
model, including turbulent pressure, was used to 
construct the mean envelope model. Frequencies 
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Fig. 4 Frequency residuals between solar models calcu- 
lated with a nonlocal generalization and the 'standard' lo- 
cal mixing-length formulation. The solid curve (NL.a-L.a) 
shows the adiabatic radial oscillation frequency shifts, un- 
der the assumption 5p t = 0, caused by the turbulent pres- 
sure pt in the mean envelope. The dashed curve (NL.na- 
NL.a) is the frequency shift caused by nonadiabaticity and 
effects of including consistently 5p t . The overall frequency 
shift (NL.na-L.a) is plotted by the dot-dashed curve (from 
Houdek 1996). 

were computed in the adiabatic approximation 
assuming 5p t = 0. 
NL.na The mean envelope model was constructed as in 
NL.a. Nonadiabatic frequencies were computed 
including consistently the Lagrangian perturba- 
tions of the convective heat flux 5F C and turbu- 
lent momentum flux 5pt- 
Additional care was necessary when frequencies between 
models with different convection treatments were com- 
pared, such as in the models L.A and NL.a. In order to iso- 
late the effect of the near-surface structures on the oscilla- 
tion frequencies the models had to posses the same stratifi- 
cation in their deep interiors. This was obtained by requiring 
the models to lie on the same adiabat near the base of the 
(surface) convection zone and to have the same convection- 
zone depth. Varying the mixing-length parameter a = H/i 
(I is the mixing length) and hydrogen abundance by iter- 
ation in model L.a, the same values for temperature and 
pressure were found at the base of the convection zone than 
those in models NL.a and NL.na. The radiative interior of 
the nonlocal models NL.a and NL.na were then replaced by 
the solution of the local model L.a, and the convection-zone 
depth was calibrated to 0.287 R© (Christensen-Dalsgaard, 
Gough & Thompson 1991). Further details of the adopted 
physics in the model calculations can be found in Houdek 
etal. (1999). 

The outcome of these calculations is shown in Fig. [4] 
As for the hydrodynamical simulations (Fig.[2]l the effect of 
the Reynolds stresses in the mean structure decreases the 
adiabatic frequencies (NL.a-L.a, solid curve) for frequen- 
cies larger than about 2 mHz, though the maximum deficit 
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Fig. 5 Frequency residuals between models for the solar- 
like star ryBoo. The model calculations and line styles are 
as in Fig.g](from Houdek 1996). 



of about 12 /zHz is smaller than in the hydrodynamical sim- 
ulations. The effects of nonadiabaticity and Sp t (NL.na-L.a, 
dashed curve), however, lead to an increase of the mode 
frequencies by as much as ~ 9 /iHz, nearly cancelling the 
downshifts from the effect of p t in the mean structure, as il- 
lustrated by the dot-dashed curve (NL.na-NL.a). If the pos- 
itive frequency shifts between models NL.na and L.a (dot- 
ted curve) are interpreted as the nonadiabatic and momen- 
tum flux corrections to the oscillation frequencies then their 
effects are to bring the frequency residuals of the hydro- 
dynamical simulations (Fig. [2]) in better agreement with the 
data plotted in Fig.Q] The effects of the near-surface regions 
in the Sun were also considered by Rosenthal et al. (1999) 
and Li et al. (2002) based on hydrodynamical simulations. 

A similar conclusion as in the solar case was also found 
for the solar-like star r\ Boo by Christensen-Dalsgaard et al. 
(1995), Houdek (1996), demonstrated in Fig. [5] and more 
recently by Straka et al. (2006). 

The near-surface frequency corrections also affect the 
determination of the modelled mean large frequency sepa- 
ration Av := (v n +u — v n i) (angular brackets indicate an 
average over n and I). In both models for the Sun and for 
77 Boo the resulting corrections to Av are about -1 /iHz. Al- 
though this correction is less than 1 % it does affect the de- 
termination of the stellar radii and ages from the observed 
values of Av and small frequency separation 5v$2 in distant 
stars. A simple procedure for estimating the near-surface 
frequency corrections was suggested recently by Kjeldsen 
et al. (2008), based on the ansatz that the frequency shifts 
can be scaled as a(v/vo) b (Christensen-Dalsgaard & Gough 
1980), where vq is a suitable reference frequency, a is de- 
termined from fitting this expression to the observed fre- 
quencies, and b is obtained from the solar data (see also 
Christensen-Dalsgaard, these proceedings). It is to be hoped 
that the analyses of high-quality data for a broad range of 
stars, observed from the space mission Kepler, will con- 
tribute to a better understanding of the near-surface effects. 
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Fig. 6 Radial pulsation mode properties of solar-type os- 
cillations for a model of the red-giant star £ Hydrae. The 
computations adopted Gough's (1977) time-dependent con- 
vection model. Top: theoretical mode lifetimes as a function 
of pulsation frequency. The horizontally dashed line indi- 
cates a mode lifetime of 15 days. Bottom: modelled oscilla- 
tion mode heights as a function of pulsation frequency (from 
Houdek & Gough 2002). 



3 Mode lifetimes in red giants 

The first convincing detection of solar-type oscillations 
in a red-giant star was announced by an international 
team of astronomers (Frandsen et al. 2002) for the star 
£ Hydrae (G7III) with velocity oscillation amplitudes of 
about 2ms~ 1 . Houdek & Gough (2002) calculated mode 
properties for the red-giant star £ Hydrae and reported 
velocity amplitudes that were in good agreement with 
the observations. Moreover, using Gough's (1977ab) time- 
dependent convection model, these authors estimated the- 
oretical mode lifetimes and reported for the most promi- 
nent p modes a lifetime r of about 15-17 days (see Fig . [6j> . 
This prediction was, however, later challenged by Stello et 
al. (2004, 2006), who developed a new method for mea- 
suring mode lifetimes from various properties of the ob- 
served oscillation power spectrum and reported a measured 
mode lifetime of only 2-3 days for the star £ Hydrae. This 
is in stark contrast to the predicted values by Houdek & 
Gough (2002), a discrepancy that needs to be understood. 
In the new method by Stello et al. it is assumed that the 
observed oscillation power spectrum is dominated by radial 
modes, because nonradial modes have larger inertiae / and 
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Fig. 7 Radial pulsation mode properties of solar-type os- 
cillations for a red-giant model in the core helium-burning 
phase. The computations were carried out by Dupret et 
al. (2009) using the convection model by Grigahcene et 
al. (2004). Top: theoretical mode lifetimes as a function 
of pulsation frequency. The radial mode lifetimes decrease 
monotonically with frequency. Bottom: modelled oscilla- 
tion mode heights as a function of pulsation frequency 
(adopted from Dupret et al. 2009). 



consequently lower amplitudes V (e.g. Dziembowski et al., 
2001, Christensen-Dalsgaard 2004). The maximum peaks 
in the observed power spectrum, however, which are actu- 
ally mode heights H and related to the velocity amplitudes 
V via the relation V 2 = rjH/2 (e.g., Chaplin et al. 2005), 
in which r\ = t _1 is the damping rate in units of angular 
frequency, are independent of the pulsation mode inertiae. 
Consequently the pulsation heights H of nonradial modes 
could have similar values to those of radial modes of similar 
frequencies, provided the modes are resolved (e.g. Dupret et 
al. 2009). A first hint, shedding some light on the pulsation 
mode lifetimes in red giants, has been provided by oscilla- 
tion data obtained by CoRoT (e.g., De Ridder et al. 2009) 
from several red-clump stars, and by Kepler (Bedding et al. 
2010). These data support the possibility of the presence of 
nonradial pulsation modes (see also Hekker, these proceed- 
ings), bringing observed mode lifetimes and theoretical pre- 
dictions in better agreement (see also Carrier et al. 2010). 

Fig.|6]illustrates the predictions of radial mode lifetimes and 
oscillation heights for a model of the red giant star £ Hydrae. 
As in the Sun the predicted mode lifetimes in £ Hydrae show 
a pronounced plateau about the maximum pulsation height 
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near 105 //Hz, a frequency value that is in good agreement 
with the observations by Frandsen et al. (2002). 
Red-giant stars may display different frequency patterns in 
the oscillation power spectrum, from regularly to very com- 
plicated patterns, depending on the density contrast between 
the core and the envelope (e.g., Dupret et al. 2009). Evo- 
lutionary calculations suggest that £ Hydrae is most likely 
in the core helium-burning phase (Teixeira et al. 2003; 
Christensen-Dalsgaard 2004). Dupret et al. (2009) com- 
puted mode lifetimes and pulsation heights for several red- 
giant models in different evolutionary phases. Mode prop- 
erties of radial pulsations for a red-giant model in the core 
helium-burning phase, calculated by Dupret et al., are dis- 
played in Fig. [7] Although the results in Figs[6] and [7j are 
for different models, though both models are in the core 
helium-burning phase, it is here assumed that their pulsation 
properties can be compared qualitatively. The most interest- 
ing qualitative difference between Figs[6] and |7j is the fre- 
quency dependence of the radial mode lifetime (top panels) 
showing a monotonic decrease with frequency in the results 
by Dupret et al. This could be a result of having adopted 
different convection formulations in the model computa- 
tions. Also, the frequency dependence of the pulsation 
mode heights (bottom panels) are different. In the stability 
computations by Dupret et al. (Fig.|7j the turbulent fluxes 
were modelled according to the time-dependent convection 
formulation by Grigahcene et al. (2004), which is based 
on Unno's (1967) formulation. Houdek & Gough (2002) 
adopted Gough's (1977ab) time-dependent nonlocal con- 
vection formulation in the model computations for £ Hydrae 
(Fig. [6}. Stability computations of red-giant oscillations 
were also recently addressed by Xiong & Deng (2007) us- 
ing Xiong's (1989) time-dependent convection formulation. 
Some of the model differences between Unno's (1967) and 
Gough's (1977a) convection formulation were discussed by 
Houdek (2008). It is to be hoped that the high-quality data 
from, for example Kepler, will not only help modellers to 
calibrate but also to improve convection models for pulsat- 
ing stars. 
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